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A search for time-integrated CP violation in D° h {h = K, tt) decays is presented using 
0.62 fb"^ of data collected by LHCb in 2011. The flavor of the charm meson is determined by the 
charge of the slow pion in the Z)*+ — > D'V^ and D*^ — > D^\^ decay chains. The difference in 
CP asymmetry between D" K+ and D" tt'tt+, AAcp = Acp{K' K+) - Acp(7r"7r+), 

is measured to be [-0.82 ± 0.21(stat.) ± 0.11(syst.)] %. This differs from the hypothesis of CP 
conservation by 3.5 standard deviations. 
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The charm sector is a promising placc^ to probe for 
the effects of physics beyond the Standard Alodel (SM). 
There has been a resurgence of interest in the past few 
years since evidence for mixing was first seen [f, 2]. 
Mixing is now well-estabhshcd [3] at a level which is con- 
sistent with, but at the upper end of, SM expectations [4]. 
By contrast, no evidence for CP violation in charm de- 
cays has yet been found. 

The time-dependent CP asymmetry Acp{f; t) for D'^ 
decays to a CP eigenstate / (with f = J) is defined as 



AcpU; t) 



V{D^{t)^f)-T{D^{t)^f) 



(1) 



where F is the decay rate for the process indicated. In 
general Acp{f] t) depends on /. For / = K~ and 
/ = ■7r~7r"'", Acp{f; t) can be expressed in terms of two 
contributions: a direct component associated with CP 
violation in the decay amplitudes, and an indirect com- 
ponent associated with CP violation in the mixing or in 
the interference between mixing and decay. In the limit 
of U-spin symmetry, the direct component is equal in 
magnitude and opposite in sign for and tt^tt"*", 

though the size of U-spin breaking effects remains to be 
qiiantificid precisely [5] . The magnitudes of CP asymme- 
tries in decays to these final states are expected to be 
small in the SM [5 8], with predictions of up to 0(10^"^). 
However, beyond the SM the rate of CP violation could 
be enhanced [5, 9]. 

The asymmetry Acp{f; t) may be written to first order 
as [10, 11] 



Acp{f; t) = a^jpif) + Jagd, 



(2) 



where ac'p(/) direct CP asymmetry, r is the 

lifetime, and aJJJp is the indirect CP asymmetry. 
To a good approximation this latter quantity is uni- 
versal [5, 12]. The time-integrated asymmetry mea- 
sured by an experiment, Acp{,f), depends upon the time- 
acceptance of that experiment. It can be written as 



Acp{f) = afp{f) + ^a^^, 



(3) 



where {t) is the average decay time in the reconstructed 
sample. Denoting by A the differences between quanti- 
ties for £)° — >• K~K^ and — >• 7r~7r+ it is then possible 
to write 



AAcP = Acp{K-K+) - Acp{n"n+) 
= [a^j.{K-K+) - a^i;(7r-7r+)] 



(4) 



A(t) 



ind 

-a^p. 



In the limit that A(t) vanishes, /S.Acp is equal to the 
difference in the direct CP asymmetry between the two 
decays. However, if the time-acceptance is different for 
the K~ and 7r~7r+ final states, a contribution from 
indirect CP violation remains. 



The most precise measurements to date of the time- 
integrated CP asymmetries in Z?" — > and Z?" 
7r^7r+ were made by the CDF, BaBar, and Belle col- 
laborations [10, 13, 14]. The Heavy Flavor Averaging 
Group (HFAG) has combined time-integrated and time- 
dependent measurements of CP asymmetries, taking ac- 
count of the different decay time acceptances, to ob- 
tain world average values for the indirect CP asymme- 
try of ag^ = (-0.03 ± 0.23)% and the difference in di- 
rect CP asymmetry between the final states of Aa^p = 
(-0.42 ±0.27)% [3]. 

In this Letter, we present a measurement of the differ- 
ence in time-integrated CP asymmetries between — >■ 
K-K+ and £>° tt-tt+, performed with 0.62 fb"^ of 
data collected at LHCb between March and June 2011. 
The flavor of the initial state (Z)° or D^) is tagged by 
requiring a — >■ D^n^ decay, with the flavor deter- 
mined by the charge of the slow pion (tt^). The inclusion 
of charge-conjugate modes is implied throughout, except 
in the definition of asymmetries. 

The raw asymmetry for tagged D° decays to a final 
state / is given by Araw(/), defined as 



^raw(/) 



N{D*+ £'0(/)7r+) + N{D*- £'0(/)7r7) ' 

(5) 

where N{X) refers to the number of reconstructed events 
of decay X after background subtraction. 

To first order the raw asymmetries may be written as 
a sum of four components, due to physics and detector 

effects: 

^raw(/) = Acp{f) + Auif) + Aj,{Tr+) + Ap{D*+). (6) 

Here, ^d(/) is the asymmetry in selecting the 1)° de- 
cay into the final state /, Ai){it^) is the asymmetry in 
selecting the slow pion from the D*~^ decay chain, and 
Ap{D*~^) is the production asymmetry for mesons. 
The asymmetries Aj) and are defined in the same 
fashion as Aj-aw The first-order expansion is valid since 
the individual asymmetries arc small. 

For a two-body decay of a spin-0 particle to a self- 
conjiigate final state there can be no D'^ detection asym- 
metry, i.e. Aj:,{K~ K~^) = AD(7r~7r"'") = 0. Moreover, 
^d(''''^) and Ap{D*^) are independent of / and thus in 
the first-order expansion of equation 5 those terms cancel 
in the difference ^raw 

{K~K+) - Araw(7r"7r+), resulting 

in 



AAcp = A,^^{K-K+) - ^aw(7r-7r+). 



(7) 



To minimize second-order effects that are related to the 
slightly different kinematic properties of the two decay 
modes and that do not cancel in AAcp, the analysis is 
performed in bins of the relevant kinematic variables, as 
discussed later. 

The LHCb detector is a forward spectrometer covering 
the pseudorapidity range 2 < 77 < 5, and is described in 
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detail in Ref. [15]. The Ring Imaging Cherenkov (RICH) 
detectors are of particular importance to this analysis, 
providing kaon-pion discrimination for the full range of 
track momenta used. The nominal downstream beam 
direction is aligned with the +z axis, and the field direc- 
tion in the LHCb dipole is such that charged particles are 
deflected in the horizontal (xz) plane. The field polar- 
ity was changed several times during data taking: about 
60% of the data were taken with the down polarity and 
40% with the other. 

Selections are applied to provide samples of ^■ 
Z?"7r+ candidates, with D° K^K^ or 7r^7r+. Events 
are required to pass both hardware and software trigger 
levels. A loose I?° selection is applied in the final state of 
the software trigger, and in the offline analysis only can- 
didates that are accepted by this trigger algorithm are 
considered. Both the trigger and offline selections impose 
a variety of requirements on kinematics and decay time to 
isolate the decays of interest, including requirements on 
the track fit quality, on the and vertex fit qual- 
ity, on the transverse momentum (pT > 2 GeV/c) and 
decay time {ct > 100 |xm) of the 13° candidate, on the 
angle between the D'^ momentum in the lab frame and its 
daughter momenta in the rest frame (| cos 6*1 < 0.9), 
that the trajectory points back to a primary vertex, 
and that the 13° daughter tracks do not. In addition, 
the offline analysis exploits the capabilities of the RICH 
system to distinguish between pions and kaons when re- 
constructing the 13° meson, with no tracks appearing as 
both pion and kaon candidates. 

A fiducial region is implemented by imposing the re- 
quirement that the slow pion lies within the central part 
of the detector acceptance. This is necessary because the 
magnetic field bends pions of one charge to the left and 
those of the other charge to the right. For soft tracks at 
large angles in the xz plane this implies that one charge is 
much more likely to remain within the 300 mrad horizon- 
tal detector acceptance, thus making A-£,{tt^) large. Al- 
though this asymmetry is formally independent of the Z3° 
decay mode, it breaks the assumption that the raw asym- 
metries are small and it carries a risk of second-order sys- 
tematic effects if the ratio of efficiencies of 13° — > K^K^ 
and D° — >■ 7r~7r+ varies in the affected region. The fidu- 
cial requirements therefore exclude edge regions in the 
slow pion {px,p) plane. Similarly, a small region of phase 
space in which one charge of slow pion is more likely to 
be swept into the beampipe region in the downstream 
tracking stations, and hence has reduced efficiency, is 
also excluded. After the implementation of the fiducial 
requirements about 70% of the events are retained. 

The invariant mass spectra of selected K^K^ and 
TT^TT"*" pairs are shown in Fig. 1. The half-width at 
half-maximum of the signal lineshape is 8.6 MeV/c^ for 
K'K+ and 11.2 MeV/c^ for 7r-7r+, where the differ- 
ence is due to the kinematics of the decays and has 
no relevance for the subsequent analysis. The mass 
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FIG. 1. Fits to the (a) m{K'K+) and (b) m{Tv'Tv+) spec- 
tra of Z)*+ candidates passing the selection and satisfying 
< Sm < 15 MeV/c^. The dashed line corresponds to the 
background component in the fit, and the vertical lines indi- 
cate the signal window of 1844-1884 MeV/c^. 

difference (Sm) spectra of selected candidates, where 
Sm = m(h~ h'^TT^) — m{h~h^) — m{'K^) for ft, = if, tt, 
are shown in Fig. 2. Candidates are required to lie in- 
side a wide 5m window of 0-15 MeV/c^, and in Fig. 2 
and for all subsequent results candidates are in addi- 
tion required to lie in a mass signal window of 1844- 
1884 MeV/c^. The D*^ signal yields are approximately 
1.44 X 10^ in the R- K+ sample, and 0.38 x 10^ in the 
7r^7r+ sample. Charm from 6-hadron decays is strongly 
suppressed by the requirement that the 13° originate 
from a primary vertex, and accounts for only 3% of 
the total yield. Of the events that contain at least one 
Z3*+ candidate, 12% contain more than one candidate; 
this is expected due to background soft pions from the 
primary vertex and all candidates are accepted. The 
background-subtracted average decay time of 13° can- 
didates passing the selection is measured for each final 
state, and the fractional difference A(t)/r is obtained. 
Systematic uncertainties on this quantity are assigned 
for the uncertainty on the world average 23° lifetime 
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FIG. 2. Fits to the 5m spectra, where the is reconstructed 
in the final states (a) K~ K'^ and (b) 7r~7r+, with mass ly- 
ing in the window of 1844-1884 MeV/c^. The dashed line 
corresponds to the background component in the fit. 



T (0.04%), charm from 6-hadron decays (0.18%), and 
the background-subtraction procedure (0.04%). Com- 
bining the systematic uncertainties in quadrature, we 
obtain A(t)/r = [9.83 ± 0.22(stat.) ± 0.19(syst.)] %. 
The 7r^7r+ and K^K'^ average decay time is (t) — 
(0.8539 ± 0.0005) ps, where the error is statistical only. 

Fits are performed on the samples in order to deter- 
mine A,-ii„{K^ K'^) and ^iaw(7r^7r+). The production 
and detection asymmetries can vary with and pseu- 
dorapidity 77, and so can the detection efficiency of the 
two different decays, in particular through the effects 
of the particle identification requirements. The analy- 
sis is performed in 54 kinematic bins defined by the 
and 77 of the candidates, the momentum of the slow 
pion, and the sign of px of the slow pion at the Z?*+ 
vertex. The events are further partitioned in two ways. 
First, the data are divided between the two dipole mag- 
net polarities. Second, the first 60% of data are processed 
separately from the remainder, with the division aligned 
with a break in data taking due to an LHC technical stop. 



In total, 216 statistically independent measurements are 
considered for each decay mode. 

In each bin, one-dimcnsional unbinned maximum like- 
lihood fits to the Sm spectra are performed. The signal 
is described as the sum of two Gaussian functions with 
a common mean fi but different widths ai, convolved 
with a function B(Sm;s) — Q{Sm) 6m'^ taking account 
of the asymmetric shape of the measured Sm distribu- 
tion. Here, s ~ —0.975 is a shape parameter fixed to the 
value determined from the global fits shown in Fig. 2, Q 
is the Heaviside step function, and the convolution runs 
over dm. The background is described by an empirical 
function of the form 1 — g-('5'"-'5™o)/a^ where Sm^ and a 
are free parameters describing the threshold and shape of 
the function, respectively. The D*^ and samples in 
a given bin are fitted simultaneously and share all shape 
parameters, except for a charge-dependent offset in the 
central value /i and an overall scale factor in the mass 
resolution. The raw asymmetry in the signal yields is 
extracted directly from this simultaneous fit. No fit pa- 
rameters are shared between the 216 subsamples of data, 
nor between the K~K~^ and 7r~7r+ final states. 

The fits do not distinguish between the signal and 
backgrounds that peak in Sm. Such backgrounds can 
arise from £)*+ decays in which the correct slow pion is 
found but the is partially mis-reconstructcd. These 
backgrounds are suppressed by the use of tight particle 
identification requirements and a narrow D'^ mass win- 
dow. From studies of the mass sidebands (1820-1840 
and 1890-1910 MeV/c^), this contamination is found to 
be approximately 1% of the signal yield and to have small 
raw asymmetry (consistent with zero asymmetry differ- 
ence between the and 7r^7r+ final states). Its 
effect on the measurement is estimated in an ensemble 
of simulated experiments and found to be negligible; a 
systematic uncertainty is assigned below based on the 
statistical precision of the estimate. 

A value of AA^p is determined in each measure- 
ment bin as the difference between A-ca,vi{K~ K^) and 
^raw(''r~7r+). Testing these 216 measurements for mutual 
consistency, we obtain x'^/ndf = 211/215 (x^ probability 
of 56%). A weighted average is performed to yield the 
result AAcp = (-0.82 ± 0.21)%, where the uncertainty 
is statistical only. 

Numerous robustness checks are made. The value of 
AAcp is studied as a function of the time at which the 
data were taken (Fig. 3) and found to be consistent with 
a constant value (x^ probability of 57%). The mea- 
surement is repeated with progressively more restrictive 
RICH particle identification requirements, finding values 
of (-0.88 ± 0.26)% and (-1.03 ± 0.31)%; both of these 
values are consistent with the baseline result when cor- 
relations are taken into account. Table I lists lS.Acp for 
eight disjoint subsamples of data split according to mag- 
net polarity, the sign of of the slow pion, and whether 
the data were taken before or after the technical stop. 
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FIG. 3. Time-dependence of the measurement. The data are 
divided into 19 disjoint, contiguous, time-ordered blocks and 
the value of AAcp measured in each block. The horizontal 
red dashed line shows the result for the combined sample. 
The vertical dashed line indicates the technical stop referred 
to in Table I. 



TABLE I. Values of AAcp measured in subsamples of the 
data, and the x^/ndf and corresponding probabilities for 
internal consistency among the 27 bins in each subsample. 
The data are divided before and after a technical stop (TS), 
by magnet polarity (up, down), and by the sign of for 
the slow pion (left, right). The consistency among the eight 
subsamples is xVndf = 6.8/7 (45%). 



Subsample 



AAcp [%] 



Pre-TS, up, left 
Pre-TS, up, right 
Pre-TS, down, left 
Pre-TS, down, right 
Post-TS, up, left 
Post-TS, up, right 
Post-TS, down, left 
Post-TS, down, right 



-1.22 ±0.59 
-1.43 ±0.59 
-0.59 ±0.52 
-0.51 ±0.52 
-0.79 ±0.90 
±0.42 ±0.93 
-0.24 ±0.56 
-1.59 ±0.57 



All data 



-0.82 ±0.21 



xVndf 



13/26 (98%) 
27/26 (39%) 
19/26 (84%) 
29/26 (30%) 
26/26 (44%) 
21/26 (77%) 
34/26 (15%) 
35/26 (12%) 



211/215 (56%) 



TABLE IL Summary of absolute systematic uncertainties for 
AAcp. 



The probability for consistency among the subsam- 
ples is 45%. The significances of the differences between 
data taken before and after the technical stop, between 
the magnet polarities, and between Px > and px < 
are 0.4, 0.6, and 0.7 standard deviations, respectively. 
Other checks include applying electron and muon vetoes 
to the slow pion and to the daughters, use of different 
kinematic binnings, validation of the size of the statisti- 
cal uncertainties with Monte Carlo pseudo-experiments, 
tightening of kinematic requirements, testing for varia- 
tion of the result with the multiplicity of tracks and of 
primary vertices in the event, use of other signal and 
background parameterizations in the fit, and imposing a 
full set of common shape parameters between D*"*" and 
D*^ candidates. Potential biases due to the inclusive 
hardware trigger selection are investigated with the sub- 
sample of data in which one of the signal final-state tracks 
is directly responsible for the hardware trigger decision. 
In all cases good stability is observed. For several of these 
checks, a reduced number of kinematic bins are used for 
simplicity. No systematic dependence of AAcp is ob- 
served with respect to the kinematic variables. 

Systematic uncertainties are assigned by: loosening the 
fiducial requirement on the slow pion; assessing the effect 
of potential peaking backgrounds in Monte Carlo pseudo- 
experiments; repeating the analysis with the asymmetry 
extracted through sideband subtraction in Sm instead of 
a fit; removing all candidates but one (chosen at random) 
in events with multiple candidates; and comparing with 
the result obtained without kinematic binning. In each 
case the full value of the change in result is taken as the 
systematic uncertainty. These uncertainties are listed in 
Table II. The sum in quadrature is 0.11%. Combin- 



Source 


Uncertainty 


Fiducial requirement 


0.01% 


Peaking background asymmetry 


0.04% 


Fit procedure 


0.08% 


Multiple candidates 


0.06% 


Kinematic binning 


0.02% 


Total 


0.11% 



ing statistical and systematic uncertainties in quadra- 
ture, this result is consistent at the la level with the 
current HFAG world average [3]. 

In conclusion, the time-integrated difference in CP 
asymmetry between — )■ K~K~^ and 13° — 7r~7r"'" de- 
cays has been measured to be 

AAcp = [-0.82 ± 0.21(stat.) ± 0.11(syst.)] % 

with 0.62 fb^^ of 2011 data. Given the dependence 
of AAcp on the direct and indirect CP asymmetries, 
shown in Eq. (4), and the measured value A(t)/T = 
[9.83 ± 0.22(stat.) ± 0.19(syst.)] %, the contribution from 
indirect CP violation is suppressed and AAcp is primar- 
ily sensitive to direct CP violation. Dividing the central 
value by the sum in quadrature of the statistical and sys- 
tematic uncertainties, the significance of the measured 
deviation from zero is 3.5(7. This is the first evidence for 
CP violation in the charm sector. To establish whether 
this result is consistent with the SM will require the anal- 
ysis of more data, as well as improved theoretical under- 
standing. 
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